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Abstract 

We study the softening of the edge modes of the compact ferromagnetic 
V = \ Quantum Hall edge as the strength of the confining potential 
decreases. To obtain a realistic potential we use a displaced background 
charge plane. We find a phase transition into a spin textured edge for 

lbB 



am^ < 0.012. 



1 Introduction 

Since the discovery of the Quantum Hall (QH) effect, work by Halperin j^] 
and, in particular. Wen has lead us to understand that there is novel low 
energy physics at the edge of an incompressible QH system. The predicted 
Luttinger liquid correlations at the edges have been studied extensively, 
leading to predictions that have considerable support in experiments 

In this paper we investigate the edge of the compact ferromagnetic u = 1 
QH system, when the potential confining the electrons to the sample models 
a realistic potential. This edge has two types of collective excitations bound 
to it: the much studied edge magnetoplasmon (EMP) ||2|, and the edge 
spin wave (ESW) When the confining potential softens (diminishes in 
strength), both these excitations soften and eventually their energies become 
negative, indicating that the edge reconstructs by spin texture (TEX) or 
charge density wave (CDW) formation |, |, |7|, |. The spin textured edge 
has the same spin structure cLS Si QH skyrmion M. 

An extensively used model for the confining potential is the u;-model 



1 10 1, where the potential is obtained from a linearly decreasing, in-plane 
positive background charge distribution. Unfortunately, however, the w- 
model is non-generic, and can not be used to determine, e.g., where, in 
parameter space, the phase transition into the TEX edge is. A natural 
attempt to obtain a more realistic description is to use the potential from 
a charge distribution displaced perpendicularly from the electron gas [pi]]. 
However, this gives a potential that converges to its bulk value very slowly, 
which leads to numerical problems. Here we handle this by modifying the 
potential slightly. 

We determine where the energies of the EMP and ESW modes become 
negative as functions of the strength of the confining potential and of ^ = 



1 



- the ratio between the Zeeman and Coulomb energies. This gives a 
phase diagram for the v = 1 QH edge. In particular, we find that as the 
potential softens the transition from the compact ferromagnetic v = 1 edge 
is into a TEX edge for g <gc = 0.012. 



2 The system and the excitations 

The QH system that we consider is a two dimensional electron gas (2DEG) 
with a large magnetic field B applied perpendicular to the plane of the 
electrons. We restrict the electrons to the lowest Landau level (LLL). For 
the vector potential we use Landau gauge; A = Bxy. The system has the 
geometry of a bar, with periodic boundary conditions in the y-direction. The 
circumference of this effective cylinder is L, and the number of electrons in 
the gas is N . The electron wave functions in the LLL are (with the magnetic 
length / = y/hc/eB = 1) 



1 



(1) 



The quantum numbers for the electrons are the momentum in the y-direction 
where n is an integer (this determines the position in the x- 



k = fn, 

direction) and the spin a =|, [. We take k = at the right edge; the left 
edge is then at —kp, where kp = ^{N — 1) is the Fermi momentum. The 
geometry is displayed in Fig. ||. 
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Figure 1: Geometry of our QH system. 



The Hamiltonian for our system is 

^iP^P')4+p',a('k+p~p',a'Ck+pyCk,a- ^ V {p - k , 0)cl^Cka 

(2) 



kpp'aa' ka,—kp<p<0 



ka 
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The first term describes the Coulomb interaction between electrons in the 
gas, the second is the interaction of the electrons with an in-plane com- 
pensating positive background charge. 6V{k) is an additional potential (to 
be discussed below), and the last term is the Zeeman energy. (A constant 
term corresponding to the self-interaction of the background charge density 
has been dropped.) In the LLL, the cyclotron energy ^hu)c is common for 
all electrons, and is ignored. The second and the third term together form 
the potential confining the electrons to the sample. The case 6V{k) = is 
called the ideal edge. c|,^, c^o- are the creation and annihilation operators 
for electrons with quantum numbers k, a. V{p,p') are the matrix elements 
for the unscreened Coulomb interaction U{r) = e^/r. Using (0), we find 

V{p,p') = ^,/Iexp-[^+^l dyKo{\p'y\)e-^y'+"^P-P">y'>/'. (3) 

L \ TT J-oo 

For the ideal edge, the groundstate of the 2DEG is the compact ferro- 
magnetic 1/ = 1 state 

1^=1)= n 4tIo)' (4) 

-kF<k<0 

i.e. all the spin up states within —kp < k < are occupied. 

We consider particle- hole (ph) excitations of (Q) within the LLL. They 
are characterized by two conserved quantum numbers: the excitation mo- 
mentum q and the spin s = 0, 1. The states are 

\q,s = 0)= Y1 '/'fcTcI+gT'^^Tl^ = 1)' (5) 

-g<fc<0 

\q,s = l) = '/'fci4+giCfctk = !)• (6) 

-kp<k<0 

The range of the momentum for the s = excitation is restricted by the 
Pauli principle. The wave functions (pka, and the corresponding eigenvalues 
(which determine the excitation energies) are calculated by numerical diag- 
onalisation of the Hamiltonian (^) in the ph-subspace {cJ.^Cp| = !)}• For 
s = 0, one finds the EMP mode, and for s = 1, one finds bulk spin waves as 
well as spin waves bound to the edge (ESW). 

3 The edge potential 

The electrons in the gas are confined by the edge potential. The Coulomb 
repulsion acts to expand the gas, and it is clear that if the potential softens 
(compared to the ideal edge potential), the gas will eventually undergo a 
phase transition to a new groundstate. 5V{k) is used to vary the confin- 
ing potential, with —k-^6V{k) as the corresponding force determining the 
softness. 
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A standard model for the confinement is the i/;- model [ |10[ . Here the 
confining potential is obtained from a compensating positive background 
charge density that decreases linearly from its bulk value to zero over a 
width w centered at fc = 0. w = corresponds to the ideal edge. For large 
enough w, phase transitions are observed. In this model, 5V{k) becomes 

5VM = J2^{p-k,0)pM, (7) 

p 

where 

{p/w + 1/2 -w/2<p<0 
p/w -1/2 0<p<w/2 
\p\ > w/2. 

For a typical 6Vw{k), see Fig. §a. The w-model turns out to be a non-generic 
choice for the confining potential, due to the property (514) (—oo) = ^^wiO) 

i- 

A more realistic model for the confinement may be obtained by displacing 
the w = charge distribution a distance d from the plane of the 2DEG - the 
displaced charge distribution (DCD) model. This models the displacement 
of the doped semiconductor layer present in real samples. In this case, one 



obtains the potential [11| 



1 / I,^2 r X + 1 

SVDCDik) = ^ / dxe'^""'^' 2darctan -+xln ^ vrd . (8) 

27r^/7r J-oo L d J 

(To derive this expression we calculated the total potential from two semi- 
infinite charge planes, at a point in the 2DEG plane.) 5VDCD{k) converges 
to its asymptotic value in the bulk, d, very slowly. This makes the numerical 
results dependent on the left edge of the sample. To overcome this finite 
size effect, we replace the DCD potential by 

5Vth{k) = ^{l-ts.n\i{k)). (9) 

^Vthik) is very similar to 5VDCD{k), as is seen in Fig. ^3. By construction, 
it has the same asymptotic values, SVthipo) = and 5Vth{—oo) = d, as 
the DCD potential, but the convergence is exponential. The parameter d 
is used to vary the potential softness; increasing d means increasing the 
softness since the force —k-^6V{k) = ^dcosh.~'^{k)k is linear in d and acts 
in the +k direction. 



4 Numerical analysis 

As seen from the Hamiltonian, the energies characterizing the system are 
the Coulomb energy e'^ /el (here we allow for arbitrary e, this can be included 
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(a) 
5V(k) 





Figure 2: The three different potentials, (a) SVw{k), w is approximately the width of 
the linear region around k = 0. (b) SVocoik) is the dotted curve, SVth{k) is the solid 
one. The asymptotic value as fc ^ — oo is d. 



in previous expressions by replacing with e^/e), the Zeeman energy gfisB 
and the energy from the confining potential. We define g = as the 

dimensionless parameter determining the bulk physics. Note that the ESW 
mode energy is linear in g, whereas the EMP is gapless. 

For a semi-infinite system (kp — > oo), the dispersion relation e(^) of the 
ESW mode is known for g ^ 1 Q, 



-f-±(5V{-oo)-5V{^))\' + 0{q' 



(10) 



here, 5 = 0. We now see explicitly what is special about the tD-model; the 
g^-term is unaffected by 5Vw{k). With 5Vth{k) as the potential in (p^, we 
get 5Vth{-'x) - 5Vth{^) = d/2 and 




(11) 



for the energy. Here we have added an index d to emphasize the dependence 
on the confining potential. Eq. ([ll] ) can be used to test our calculations. 
For g <C 1 at a given d, our numerical values should converge to the energy 
given by ed{q) when the system size increases. 

We want to calculate energies for a half-plane, i.e. in the limit N/L oo 
(the width of the system is kp oc ^) and L ^ oo. Upper limits for the 
parameters are set by the simulation time, which grows dramatically with 
A^. In our simulations we varied N between 200 and 800 and L between 50 
and 400, with 1.5 < N/L < 6. To find the asymptotical values, we scale in 
L/N and l/L and extrapolate to zero. (When scaling in L/N we keep 1/L 
constant and vice versa.) 

The energies of the ESW and EMP modes only change appreciable with 
L/N for small q. (This is true for all d.) At g ~ 0.6 the energies differ by 
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~ 10~^ for L/N between 0.17 and 0.5 (at constant L). For larger L/N , the 
difference increases. It also grows with decreasing q. At g ~ 0.1 for instance, 
the difference is ~ 10~^ for L/N between 0.17 and 0.5. For q < 0.4, we use 
scahng in both L/N and 1/L. For 0.4 < q < 0.6, we use 1/L scahng with the 
restriction L/N < 0.5, whereas for q > 0.6, system sizes down to L/N = 2/3 
were used. 

First we test our results at g <C 1. Fig. ^ shows the ESW energy as a 
function of L/N and 1/L at d = 0.65 and q = 0.0628, compared to the q^- 
term of (pi]). In Fig. ^ the energy is seen to approach the predicted value 



Figure 3: Scaling diagrams for the ESW at d = 0.65 and q = 0.0628. (a) The energy as 
a function of L/N, with L = 100. The data points are from simulations, and the dashed 
curve is the 0(g^)-value predicted by (pT|). (b) The energy as a function of 1/L. The 
upper data has L/N — 0.33 and the lower L/N = 0.50. The line is a least-squares fit. 



nicely. (We do not expect an exact correspondece since we have a finite q.) 
In Fig. |b we see that the effect of increasing L is marginal, certainly when 
compared to the energy variation in Fig. This analysis can be repeated 
for different values of d, and the agreement with the model in ( |TT| ) is good 
at least for 0.05 < d < 0.80. We conclude that our numerical results are 
trustworthy. 

To find the critical points where the modes go soft, we first isolate the 
approximate d and q values where this happens. Fig. ^ shows typical results 
for the dispersion relations for both types of excitations for two different 
values of d, with ^ = 0, = 600 and L = 100. 

In Fig. ^ we are close to the point where the ESW mode goes soft. 
We determine the d value more exactly and in Fig. |5| we see the minimum 
energy as a function of 1/L for two different values of d, with the minimum 
occuring at q ^ 0.48 (weakly depending on L due to the quantisation of 
q). We see that the energy varies quite strongly with L, however, a linear 
least-squares fit matches the data very well. Using this line to extrapolate 
we get dESW = 0.602 for the critical value. Note that the analytic result, eq. 
([Tl|), predicts a phase transition at g = for d ~ 0.627 (if only the g^-term 
is considered), whereas we find the transition for finite q, where ed{q) does 
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Figure 4: Dispersion relations for the two excitations, sd, g = 0, N = 
The sohd curve is the ESW mode and the dotted the EMP. In (a) d — 
d = 0.650. 



600 and L = 100. 
0.600, and in (b) 



D.006 

1/L [1-'] 



Figure 5: Minimum energy for the ESW mode as a function of 1/L. The upper data 
points are at d = 0.601, the lower at d = 0.602. The lines are least-squares fits. 



not apply, and a smaller d. We can thus safely conclude that we have found 
the correct minimum. 

In Fig. d ~ dEMP, where the EMP mode goes soft. We examine the 
energy at the minimum (q ~ 1.28) for the EMP mode as a function of 1/L 
for two different values of d in Fig. |6|. We clearly see that dsMP = 0.652 is 
the best approximation. 

We also see from Fig. |^ that the ESW mode has continued to drop 
below zero as we increased d. In order to have a CDW phase, we need the 
ESW mode energies to be positive. Since the ESW mode is linear in g, the 
energies will become positive if ^ > gc, where gc is the absolute value of 
the minimal ESW energy, at d = dEMP- Fig. |7| shows the minimum ESW 
energy (at q ~ 0.75) as a function of 1/L, for d = 0.651 and 0.652. The L 
dependence is fairly strong, but again the data is practically linear in 1/L. 
We already know that dEMP is closer to 0.652 than 0.651, and, as we can see, 
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Figure 6: Minimum energy for the EMP mode as a function of 1/L. The upper data 
have d = 0.651, the lower d = 0.652. 
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Figure 7: Minimum energy for the ESW mode as a function of 1/L. The upper data 
points are at d = 0.651, the lower at d = 0.652. The lines are least-squares fits. 



Qc only changes by ~ 10 in that interval. We conclude that Qc = 0.012. 



5 Summary 

Along the curve from {Q^dEsw) to {gc,dEMp) in gd-space there is a phase 
transition into a TEX edge. Similarly, along the line d = dEMP and g > Qc 
there is a phase transition into a CDW edge. The critical values are dESW = 
0.602, dEMP = 0.652 and Qc = 0.012. 

In our calculations we assume that the groundstate of the edge is the 
compact ferromagnetic \v = 1) state, and, therefore, we cannot determine 
anything beyond the TEX and the CDW phase transitions. At larger d, 



however, a phase that is a combination of the two is expected [12|. Fig. 
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shows the phase diagram. 

For the it;-model, using the explained numerical techniques, we obtain 
the asymptotical value gc^^, = 0.008 for w = 7.03. Thus the u)-model under- 
estimates the region where the transition is into the TEX edge. 
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Figure 8: Phase diagram in gd-spacc with the different regions corresponding to different 
groundstates. Solid lines are numerically determined, dotted ones are only topologically 
known. 
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